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Abstract 

Fluctuations in electron density and tempera- 
ture coupled through Ohm's law are studied for an 
ionizable medium. The nonlinear effects are con- 
sidered in the limit of a third order quasi-linear 
treatment. Equations are derived for the ampli- 
tude of the fluctuation. Conditions under which a 
steady state can exist in the presence of the 
fluctuation are examined and effective transport 
properties are determined. A comparison is made 
to previously considered second order theory. The 
effect of third order terms indicates the possi- 
bility of fluctuations existing in regions pre- 
dicted stable by previous analysis . 

Introduction 

Under conditions appropriate to MHD gener- 
ators and MPD arc thrusters studies have shown 
that fluctuations in electron density and/or tem- 
perature can be unstable. In both of these 
analyses only the linear effect of fluctuations 
was considered leading to the usual stability 
analysis dispersion relation from which the growth 
or damping of the fluctuation can be ascertained. 
In this paper we extend these analyses into the 
nonlinear regime in order to determine the ampli- 
tude of the fluctuation and its effect upon the 
properties of interest. 

The oscillations under study primarily arise 
due to fluctuations in electron temperature. 

These temperature fluctuations can be amplified 
through increased Ohmic heating arising as a re- 
sult of increased ionization and/or decreased col- 
lision frequency. These changes in frequency 
occur as a consequence of the initial electron 
temperature perturbation. The dispersion relation 
depicting this phenomena was derived in refer- 
ence (2) by considering the time dependent elec- 
tron density and temperature equations coupled 
through a generalized Ohm's law. 

In this paper the nonlinear effects of the 
fluctuations are considered on the basis of a 
third order quasi-linear theory neglecting mode 
coupling. Even though the theory is strictly 
valid only in the neighborhood of the threshold 
of stability, the third order terms are included 
in an attempt to increase the applicability of the 
results. The analysis is carried out for those 
fluctuations which show the maximum rate of 
growth, i.e. , first become unstable. Two limiting 
cases are considered. In the first limit, the 
ionization is assumed to be frozen. In this limit 
it is shown that the system can become unstable if 
electron collision frequencies decrease suffi- 
ciently rapidly with electron temperature. For 
the MHD generator, this mode does not appear to be 
important since operating conditions and working 
fluids of interest are such that the instability 
criterion is seldom satisfied. One notable excep- 
tion may occur, however, when the gas is strongly 
preionized prior to entering the MHD channel. In 


this case, the ionization may be sufficient that 
Coulomb collisions dominate and the instability 
criterion is satisfied. This was mentioned in 
reference (3) as a possible source of the anomo- 
lous oscillations observed in these experiments 
with preionization. In reference (3) the oscilla- 
tion was referred to as the static instability 
mentioned in reference (1) and is akin to our 
frozen flow limit. In this limit, the medium is 
assumed to be infinite and homogeneous and the re- 
sults of the quasi-linear analysis are presented 
for a fluctuation propagating in the direction of 
maximum growth. 

For the MPD arc thruster, the fluctuations 
are assuned to occur in the region of high current 
concentration in the throat of the device which is 
depicted by an annular region between two concen- 
tric cylinders as shown in figures 1 and 2. Since 
damping due to radiation and heat conduction are 
least prominent for long wavelengths, the most un- 
stable waves are those with the longest wave- 
lengths. In the case where the interelectrode 
distance is small compared to the mean circumfer- 
ence of the annular region, the longest wavelength 
is obviously nearly that of purely rotational 
propagation. The analysis of the MPD arc thruster 
is therefore simplified by considering only rota- 
tional propagation. 

In the second limiting case, the ionization 
rate is assumed infinite. This case has been ex- 
tensively analyzed with regard to MHD generators 
and was reviewed in reference (4). In that study, 
terms were retained to second order in fluctuating 
quantities. In our work, third order terms and 
some terms ignored in reference (4) as being small 
in the range of interest are retained. In the 
range of interest, i.e., the region in which sta- 
bility is reached in the presence of fluctuations 
of finite amplitude, the effect of these addi- 
tional terms is minor, at least in so far as 
present experimental data is interpretable. How- 
ever, an interesting point arises in that these 
additional terms open the possibility of insta- 
bility occurring in a region where linear theory 
predicts stability. In this case, instability can 
occur if a fluctuation of sufficient initial am- 
plitude is present. Such fluctuations may be ex- 
cited in the preionization region of the MHD duct. 

For MPD arc thrusters in the purely rota- 
tional approximation the mode occurring in the in- 
finite ionization rate limit is damped. For prop- 
agation at small angles from the rotational direc- 
tion the mode can be unstable but requires large 
values of the Hall parameter for instability to 
occur. From the analysis of reference (2) , these 
values appear to be much higher than those which 
occur in the high current throat area of present 
high pressure MPD thrusters. Therefore, it is not 
expected that this mode of instability plays an 
important role in the rotating spoke phenomena 
observed in these devices. 
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Assumptions and Governing Equations 


Generalized Ohm's Law 


Assumptions 


The generalized Ohm's Law is 


In those regimes where electrothermal dis- 
turbances predominate, the following assumptions, 
simplifications, and restrictions are appropriate: 


= oE 


v * X ^0 


( 4 ) 


1. The analysis is restricted to sufficiently 
short wavelengths such that the zeroth order prop- 
erties do not vary appreciably over this distance, 
i.e., spatial gradients in zeroth order quantities 
are ignored. 

2. Only fluctuations in space and/or time of 
the electron number density and temperature are 
considered. This is reasonable since the rela- 
tively lighter and more mobile electrons respond 
to disturbances of greatly different frequencies 
than do the heavier and less mobile neutral atoms 
and ions . 

3. Only the propagation of magnetohydrody- 
namic disturbances is considered. This restric- 
tion implies neglecting the displacement current 
density in Maxwell's equations, ignoring induced 
magnetic field relative to the applied constant 
magnetic field, and assuming the plasma to be 
quasi-neutral. 

4. Ion and neutral particle flow velocities 
are equal, i.e., ion slip is ignored. Further- 
more, all heavy particle flow velocities are as- 
sumed equal to V^, the gas flow velocity. 


** * 1 — 

* = * + tt 7p e < 5 > 

e 

* 

where E = j£ + Vq x is the electric field in 
the frame of reference moving with the gas flow 
velocity Vq in the applied magnetic field Bq, 

and (l/n e e)Vp e is the force due to electron tem- 
perature and density gradients. 


The electrical conductivity, o, is given by 


a 


2 


n e 
e 

m v 
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( 6 ) 


where the total electron momentum collision fre- 
quency is 


v = v + v + v . 
ec es ei 


(7) 


The species collision frequencies are allowed the 
following general temperature dependence 


v = A n T 
ea ea a e 


( 8 ) 


5. Only propagation in planes perpendicular 
to the applied magnetic field is considered. 

6. Ion and neutral particle temperatures are 

assumed equal and equal to T , the gas tempera- 
ture. 1 


where, a = c, s, or i corresponding to electron 
collisions with carrier gas atoms, seed atoms, or 
seed ions, respectively. A ea and a ea are arbi- 
trary constants chosen to fit average collision 
frequency data over the temperature range of in- 
terest . 


For the above limitations the problem is com- 
pletely determined in terms of the fluctuations of 
electron density and temperature and their cou- 
pling through Ohm's law and Maxwell's equations. 

Maxwell's Equations 


Electron Energy Equation 


The electron energy equation is 




T 

e 


+ ^0 



• 7 


kT 
e 
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From Maxwell's equations under the assumption 
of quasi-charge neutrality and neglecting the in- 
duced magnetic field, 

7 x E = 0 (1) 


= -2v m n g § k(T e - V + i ■ E* 

-( kT i + ! kT e)(lr n e + ^0 • ’ n e) (9) 


7 • 3 - = 0 

Electron Continuity Equation 


( 2 ) 


It is assumed that in the region of interest, 
the ionization is dominated by electron-neutral 
atom ionizing collisions and by three body recom- 
bination. The electron density equation then 
takes the form 


37 n e + ^0 


7n 


n n v. 
e s i 


3 
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(3) 


It is noted that this equation depends upon the 
gas flow velocity 7 q , rather than the electron 
flow velocity, 7g. This arises as a result of 
assumption (4) and the condition of quasi-charge 
neutrality . 


where v is the collision frequency for energy 
transfer given by 

m m m 

v = — v + — v + — v . (10) 

m m ec m es m. ex 

c s i 

and the 7T g contribution to the heat conduction 
term and radiation losses are neglected. These 
terms contribute damping terms to the dispersion 
relation which vary inversely with the wavelength 
relative to the wavelength independent elastic 
collision damping. Since, in quasi-linear theory, 
one is primarily restricted to the region in the 
neighborhood of the stability threshold, only 
those waves which are most unstable are of inter- 
est, i.e., the longest wavelengths. It is there- 
fore assumed that at these wavelengths, the wave- 
length dependent terms are ignorably small and 
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hence for simplicity are not included in the anal- 
ysis. The question of the relative importance of 
these terms has been extensively investigated for 
MHD generator type plasmas in references (5) and 
and (6) . 


Ouasi-Linear Theory 

The nonlinear problem is now treated by con- 
sidering small departures from the stability 
threshold so that pertinent quantities can be 
written in the form 

f (r , t) = <f> + f(r,t) (11) 

where f is the oscillatory function for which the 
time averaged value defined by the brackets (f) is 
zero and (f) is the average value of the quantity 
f. Furthermore, in keeping with our initial 
assumption (f) is spatially independent. The re- 
striction of small departures from steady state 
then requires 


f « (f> (12) 

so that a solution may be formed as a power series 
expansion in f. In the cases considered here, the 
expansion is carried out to third order in the 
product of i's. Furthermore, since we are con- 
sidering small departures, only the fundamental 
mode is considered, i.e., mode coupling between 
harmonics or other fundamental modes is ignored. 

The perturbation function is then taken to be 
of the form 


f(r_, t) = f s (t)sin(jt • r_ - mt) 

+ f c (t)cos(^ • £ - ut) (13) 

where and m are real, 1 • Bo = 0 (only con- 
sider propagation in plane perpendicular to Bq) , 
and the long term time dependence of the perturba- 
tion function is contained in the coefficients f s 
and f c which are assumed to be slowly varying 
over a period of the oscillation, t, so that 


< f s> 5 X / t+T f s (t '> dt ’ s f s (t) | 

i r t+T 

< f c> S 7 / f c^’> dt ’ " f c (t) 
•'t J 


(14) 


The problem of interest is then one of determining 
under what circumstances, if any, that the ampli- 
tude of the fluctuation reaches a steady state 
value. 


The general quasi-linear analysis of the set 
of equations (1) to (10) is complicated by the two 
time dependent equations (3) and (9) . This allows 
two fundamental modes to propagate simultaneously. 
Therefore, in order to simplify the analysis, we 
consider two limiting cases. In the first case, 
oscillations in the electron number density are 
ignored. This is justified if either the ionizable 
material is completely ionized or the frequency of 
the oscillation is much greater than the ionization 
frequency. The electron number density like the 
gas dynamic properties is then specifiable by equa- 
tions peculiar to the problem of interest. 


In the second case we assume the ionization 
rate is infinite. In this case the left hand side 
of equation (3) can be neglected. The electron 
number density is then specified by 


2 

n 

e 




(15) 


It is further assumed that the electrons are ini- 
tially in Saha equilibrium so that we are only 
concerned with small fluctuations from equilibrium. 
Under this condition, the principle of detail bal- 
ance^) can be invoked to relate the ionization 
coefficient, Vj, to the recombination coefficient, 
v r , through the equilibrium constant. Then 



This case has been extensively analyzed with re- 
gard to MHD generators . (4) this paper, we ex- 

tend this second order analysis to third order and 
also consider it with regard to high pressure MPD 
arc thrusters. 

Frozen Ionization Rate Limit 


This is the case where electron number density 
fluctuations are ignored. Before proceeding to the 
quasi-linear analysis, we consider the results 
previously obtained(2) from linear stability anal- 
ysis. For these results extended to include a 
neutral carrier gas as well as an ionizable species 
the condition that the system be unstable is 



X 




a 


es 


V esO 

v 0 


+ a 


ei 



< 0 


(17) 


Let us first discuss the consequences of this 
stability criterion for combustion and nonequilib- 
rium MHD generators and then consider the MPD 
thruster case. In the combustion generator in- 
elastic collisions with the numerous molecular 
species in the combustion gas work fluid make it 
impossible to raise the electron temperature much 
above the gas temperature. In this case then, the 
factor (1 - T 0 /T e0 ) in the inequality (17) is in 
general, too small for the instability criterion 
to be satisfied. 


In the nonequilibrium MHD generator operating 
on inert gases seeded with alkali metals the goal 
is to obtain an elevated electron temperature. An 
electron temperature approximately twice the gas 
temperature is typical. There, even for the most 
unstable direction of propagation, i.e., 

L ’ 3-q = the unstable condition requires the 


coefficient 


ecO 


ec v n 


+ a 


esO 


+ a 


eiO 


< - 1 . 


In general, over the temperature range and for 
neutral gases of interest, the temperature expo- 
nentials a ec and a gs are positive or at least 
only slightly negative (e.g., see Maxwell-averaged 
cross sections for electron temperatures near 
3000° K presented in ref. (8)). Therefore, in 


3 



L ■ <0 ( 23 ) 
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order for the coefficient to be < -l, the plasma 
must be Coulomb collision dominated since 
a e £ - -3/2. Since under normal operating condi- 
tions of MHD generators the collision frequency is 
neutral dominated, it appears that the range of 
importance of this mode of instability for these 
devices is small and can be avoided without great 
difficulty. The one exception to this could occur 
when the gas is preionized by strong electric 
fields prior to entering the MHD channel. In this 
case, high degrees of ionization can be reached 
and the resulting instability could propagate into 
the generator region. This has been mentioned in 
reference (3) as a possible source of the anomo- 
lous fluctuations observed in these experiments. 

On the basis of the above discussion and the 
fact that MPD arc thrusters in general operate in 
the Coulomb collision dominated regime, the quasi- 
linear analysis in the limiting case of frozen 
ionization is developed only for v c v e ^. 


(}) = (a)(E*) - x £ 



It is convenient at this point to define for 
later use the effective values of the conductivity 
and Hall parameter. These quantities are defined 
by 

<*> E °eff<^>- B eff<3> x i ( 24 ) 

It is obvious from equation (24) that 


eff 


<4> 2 

’eff - ^ 

i • (<i> x <e*>) 


(25) 

(26) 


We first consider the equations governing 
the average and fluctuating electric current 
densities. Substituting quantities of the form 
of equation (11) into equation (4) , we obtain 


Substituting equation (23) into equations (25) 
and (26) 
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eff 


1 + Z 


(27) 


4 = (<0> + 5)«e*> + £**) - «B> + B) (<J-> + l) x £ 

(18) 


= .<g) +_ ZT. 
eff 1 + Z 


(28) 


where to third order (]2 ) = (E, ), £ is a unit 

vector in the direction of Bq and the Hall 
parameter, S, is defined by 



(19) 


Averaging equation (18) then gives an expression 
for the average current density. 


<j> = <°Xl*> + <!**> - <B)<j> x £ - <Bd> x £ 

( 20 ) 


The fluctuating part of the current density 
is then obtained by subtracting equation (20) from 
(18) , forming the triple cross product 

x ^E x (j^) , and using equations (1) and (5), 
i.e. , E. x I** = 0, and equation (2), i.e. , 

_E ■ J = 0 to obtain 


where 


*— =- (29) 

(1 + t ) 



Much of the above manipulation was from hind- 
sight and for the purpose of obtaining the equa- 
tions for the average current density, equa- 
tion (23), and the fluctuating current density, 
equation (21) , solely in terms of the fluctuating 
electrical conductivity, a, and the fluctuating 
Hall parameter, B. *rhe above equations are com- 
pletely general, i.e., not restricted to either of 
our limiting cases. Specialization to these lim- 
itingcases arises through the dependence of 3 
and B upon n e and T e in these limits. 



I = [oE x (E x (e*)) - BE x [ix «^) x £) ]J 

( 21 ) 

The fluctuating part of the field, E , can 
be determined from the fluctuating current density 
equation obtained by subtracting equation (20) 
from (18), forming the dot groduct E. • |, using 
equation (1) so that E = — V<J> , and equation (2), 
i.e., E_ • j = 0 to obtain 

<0>E 2 r = -0E(E • <E*>) + EBE • (<i) x b) 

+ E<B>1 • (4 X b) 

+ higher order terms (22) 

Equation (20) can now be simplified by using equa- 
tion (22) to eliminate ^oE* / and equation (21) to 
eliminate We then obtain 


For the case under consideration 


n = a = a 
e ec es 


°; « el = -3/2 


then from equations (6) , (7) , and (8) 
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<T 2 > 
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16 


( T e> 


Q * (T 2 ) 

9 e \ e' 

16 


(31) 


From equations (19) , (7) , (8) , (31) , and (32) 


( 32 ) 
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6 0 °0 


B 3 

W ‘ <34> 

where 0g and ^3 are defined by equations (6) 
and (19), respectively, evaluated at the average 
electron temperature, (T g ). 

We now consider the energy equation (9) by 
eliminating 4 • E* in terms of (l/o)4 2 by tak- 
ing the dot product of j with equation (4) . It 
should then be obvious from equations (21), (23), 
(32), and (34) that the energy equation can be ex- 
panded in a power series in T e /(T e ). Carrying 
out this expansion to third order in f e /^T e ) and 
averaging the resulting equation we obtain, after 
considerable algebraic manipulation, the following 
equation for the average electron temperature 


e0 < T «> * ' 2 IT "e0”0 2 k< < T «> ‘ 1 
a 
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L »«>1 

<0 : 

1 

8 < T e>l 

" (< T e>“ T (? 
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/x 2 N " 

+ ^T 

<4> 2 + (A • 

l 

/i» 2 9 ' 


(35) 


At this point, it is appropriate to consider 
the essence of the analysis. First, it is to de- 
termine under what conditions, if any, the ampli- 
tude of a fluctuation which is linearly unstable 
can reach a steady state. The second purpose is 
to determine the values of the measured properties 
of the system in the presence of fluctuations. 
Since most measured properties are obtained with 
instruments that cannot resolve oscillations on 
the time scale of the fluctuation, they only meas- 
ure time averaged properties (ideally) . 


It may be noted that for the third order the 
theory considered here, the pertinent equations for 
the average values, e.g. , equation (31), depend 
upon the average of the square of the amplitude, 
^2^ rather than the amplitude itself. Therefore, 
it is convenient to cast the equation for the 
amplitude, T e , in this form. The equation for T e 
is obtained by subtracting equation (35) from equa- 
tion (9), expanding to third order in T e ^ T e^’ 
multiplying by T e and averaging to obtain 


2 kn eo( T e) 2 at ^e) ' 2 m n e0 v 0 2 k( ( T e) V 


f 3 35 3 ^e) (*e) . . 

15 3 (O > 

| 2 ' 16 2 <T e > 2 ( < T e) - V 

8 2 ^f 

r~ „ “1 
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3 21 Ve> 

<^> 2 

2 32 <T e ) 2 


- 2 


2 + 51 ^,e)_ 


2 ' ^ < T e> 2 


72 |<0 

(36) 


which closes the set of equations describing the 
average properties of the system in terms of the 
average of the square of the amplitude of the 
fluctuation, ^t|^. 

Before proceeding to analyze the above equa- 
tions, another property of interest, namely, the 
frequency of the fluctuation is determined. This 
quantity is obtained by forming the equation for 
the amplitude of the fluctuation, T e , as in the 
above derivation, but instead of multiplying_by 
T e and averaging, we multiply by l/i 2 • VT e 
and average. The frequency in the frame of refer- 
ence moving with the gas flow velocity is then 

* A • <4> 

“ = “ + - ’ ^0= (37) 

which was the expression previously determined 
from linear stability analysis. ^ 2 ) in other 
words , the frequency of the oscillation depends 
upon the amplitude of the oscillation only through 

<*>• 


The Effect of Fluctuations on Ohm's Law . The 
conductivity and Hall parameter are often indi- 
rectly determined in experiments by first measur- 
ing the electric current density and electromag- 
netic fields. Then Ohm's Law is applied to deter- 
mine the conductivity and Hall parameter. In the 
presence of fluctuations it is then obvious from 
equation (24) that the measured values are the 
effective values. Clearly, a discrepancy between 
theory and experiment arises when one calculates 
the conductivity and Hall parameter in their clas- 
sical form, i.e., equations (6) and (19) evaluated 
at the measured electron temperature. It is 
therefore of interest to compare the ratio of ef- 
fective to classical values of these parameters. 


From equations (27) , (29) , (30) , and (31) we 
obtain 



In figure 3, this ratio is plotted as a function of 
the amplitude function for various 

angles of propagation. The interesting point to 
note is that the effective conductivity is en- 
hanced by the presence of the fluctuation for 
values of |i| > /5 but is reduced for values lesc. 
than /E. 


From equations (28), (29), (30), (33), and 
(31) , we obtain 
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eff 
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<Q\,9 x ( T e) 

<T e )V Sl+T 2 )^) 2 
4 a + t 2 ) <T e ) 2 


(39) 


Comparison of equations (39) and (38) shows that 
in the limiting cases x = 0 ai.d jx| = 03 > 


5 



At intermediate values of 


depends 


upon the value of Bg. Therefore, the dividing 
point between enhancement and degradation of the 
Hall parameter is not simply a function of the 
angle of propagation as it was for the conduc- 
tivity ratio. The plot of 3 e ff/3o for the lim- 
iting cases is therefore identical to the a e ff/ c o 
plot in figure 3, except that the neutral point 
Beff/00 = 1 occurs for values of x given by 




These values of x are plotted in figure 4 for 
various values of Bg. The angles of propagation 
for which B e ff = Bg are restricted to 
0< x < /5 and £ x <_ -/J. The asymptotic 
values x = occur in the limit 8g 00 . At 

the limiting point equality (40) again holds as 
can be seen by taking the limit gg -*• “ in equa- 
tion (39) and comparing the results to equa- 
tion (38) . 

The Effect of Fluctuations for the Case of 
Maximum Growth Rate . We are interested in deter- 
mining whether a steady state can exist in the 
presence of an infinitesimal oscillation, and if 
so, under what circumstance the oscillation occurs 
and what its amplitude is as a function of experi- 
mentally measured parameters. 

We first consider the amplitude equation (36) 
for the direction of propagation which results in 
the maximum growth rate of the disturbance. Since 
in this direction, the fluctuation is first un- 
stable and grows most rapidly, it is assumed that 
for small displacements from the stability thresh- 
old this is the observed disturbance. For the 
Coulomb dominated regime considered here 
(a ec = a es = 0, a e i = -3/2) this direction, as can 
be seen from the linear stability term of equa- 
tion (17) , occurs for propagation perpendicular to 
the average current density, i.e. , £ • = 0. 

For this case, the amplitude function is given by 


3 29 (z_ T 0 \ 
2 8 \29" fT)) 


where we have used equation (35) with 

3/3t (T e ) =0 to eliminate (l/(o'>) and have 

neglected terms greater than second order in 

. We note that the linear term is just 
the stability condition obtained from the linear 
analysis and obviously is damping until the in- 
stability point, Tg/^T e ^ £. 2/3 , is reached. Be- 
yond this point the linear term causes a growth in 
the amplitude of the oscillation until the non- 
linear term, which is damping for Tg/^T e ^ > 2/29, 
offsets the linear growth term. A steady state is 
reached when 


1 29 (g_ T 0 \ 

2 8 \29 ~ <T e )j 


Obviously in the case in which T g/(T e ) < 2/29 
the nonlinear term contributes to the growth and 
the thermal perturbation is unconditionally un- 
stable. However, it is obvious from equation (43) 
that as T 0 /<T e ) 2/29, <T 2 >/<T e > 2 + - which vi- 

olates the quasi-linear condition of small ampli- 
tude oscillations, i.e., ^2^/^T e ^ 2 << 1. There- 
fore, this point is not considered further. 

The amplitude function, (i|)/^T e ^ 2 , as 8 : *- veri 
by equation (43) , can be calculated by determining 
the temperature ratio T g/(T e ^ from equation (35). 
For l • = 0, the steady state condition 

3/3t^T e ^ = 0, and either the MPD thruster geometry 
shown in figure 2 or the MHD Hall generator con- 
figuration, this equation reduces to 





i + <e> 2 ) 


Hkl 


is the ratio of the ion drift energy to its random 
kinetic energy, E = |\E* and Bg are the magni- 
tude of the applied electric and magnetic fields , 
respectively. In deriving equation (44), we have 
also used equations (23), (3l) , and (33). 

Equations (43) and (44) can be solved simul- 
taneously for given values of the Hall parameter 
and the parameter R to determine the ampli- 
tude function (TeVO-e/ * These values are shown 
in figure 5, where the amplitude function has been 
plotted as a function of the Hall parameter (b^ for 
various values of the parameter R. The first 
thing to be noted from this figure is that for a 
given value of R, there is no instability, i.e., 
no physically exceptable solutions exist, until a 
critical value of is reached corresponding to 
the onset of instability. This critical value of 
<e> increases with decreasing values of R. Also 
to be noted from figure 5 is that the amplitude 
function reaches an asymptotic value as (8) ■», 

the magnitude of which depends upon the parameter 
R. 

The Effect of Fluctuations for the Case of a 
Purely Rotational Disturbance . In the previous 
section the medium was taken to be infinite and 
homogeneous. As explained in the Introduction in 
the case where the disturbance is confined to the 
annular region between two concentric cylindrical 
electrodes, the longest wavelength and hence the 
most unstable wavelength (due to radiation and heat 
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conduction damping) is approximately that which 
corresponds to pure rotation. In this section, we 
therefore consider the effect of a purely rota- 
tional fluctuation. This was the case considered 
in reference 2. 


The amplitude equation (36) to second order 
in 01^0- eV^ reduces, after considerable alge- 
braic manipulation, to 



3 29 2 T 0 23 2 53 2 T 0 3_ ft 4| 

2 8 [29 29 0 “ 29 0 29 P 0 


, g>) < r l> 

« * e o> < t «> 2 J <?$ 


(46) 


Unfortunately the dependence of equation (46) upon 
Tg/^T e } and B^ is such that it is not possible 
to determine a steady state amplitude of the fluc- 
tuation for the same assumptions used in the line- 
arized theory of reference (2) . In order to ob- 
tain the values of the local gas dynamic proper- 
ties required in the theory from the available ex- 
perimental data, a number of assumptions had to be 
made in the analysis of reference (2) . One of 
these was to take the local gas temperature, Tq, 
to be sufficiently small so that it could be ig- 
nored in the onset region where the frequency of 
the oscillation was evaluated. If this is done 
here, the nonlinear term of equation (46) becomes 


3 29 
2 8 


[2_ + 23 fi 2 3_ ft 4l 
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, <o 

<1 * «J> 


(47) 


which is positive and hence destabilizing even for 
the maximum value of Bq for which the linear 
term in equation (46) is destabilizing, i.e., 

Bq = 2. Therefore, a more complete steady state 
solution is required and/or local gas properties 
must be measured experimentally in order to ascer- 
tain the validity of this assumption. 


However, a second interpretation is possible. 
Since, in reference (2), the rotational frequency 
was calculated at what was assumed to be the up- 
stream onset point, it could well be that the 
amplitude of the oscillation is growing at this 
point, reaching its quasi-steady value downstream 
of the point in question. In this case, the 
amplitude may indeed be small at the point at 
which the frequency was calculated so that non- 
linear effects can be ignored. The effect of non- 
linearity then manifests itself downstream in the 
growth of the oscillation. 


Of further interest is the fact that the 
linear and nonlinear coefficients in equation (46) , 


C ° S 2 ‘ 3 < T e> 
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respectively, can be positive or negative, depend- 
ing upon the magnitude of Tq/(T 6 ) and Bq. In 
figure 6, the curves of Cq = 0 and " " 


Cj = 0 


plotted as functions of Tq/^T^ and B^. For both 
coefficients, points (Bq, T o/(T e ^) above their re- 
spective curves result In negative values of the 
coefficient and hence damp the fluctuation. How- 
ever, for points (8$, T o/( T e )) below their respec- 
tive curves, the coefficients are positive so that 
the terms contribute to the growth of the fluctua- 
tion. As a result, the Tq/^T 6 } against Bq plane 
is divided into four regions having the following 
characteristics : 


1. Region I - Both coefficients are damping 
and hence fluctuation is unconditionally damped. 

2. Region II - Linear term damping, nonlinear 
term growing. In this case the usual linear sta- 
bility theory would predict a damped disturbance. 
However, the growth contribution of the nonlinear 
term indicates that, for disturbances with initial 
magnitude of the amplitude function 

<*iy< T e y > ICq/C-J , the disturbances grow un- 
stably. Therefore, in this region fluctuations 
are damped for |^/^T e ^ < | Cq/C^ | and grow un- 
stably for <T2)/(r e ) 2 > (Cq/C-J . 

-3. Region III - Both coefficients are positive 
and hence any disturbance grows unstably. 

4. Region IV - Linear term growing, nonlinear 
term damping. Disturbance grows or damps until a 
stable equilibrium point is reached with amplitude 
function <T|>/<T e ) 2 = | Cq/C^ | . 

Whether the above regions exist under actual 
experimental conditions only further experimenta- 
tion can answer. Furthermore, the effect of higher 
order terms in those regions in which the theory 
predicts the unstable growth of the disturbance is 
not known. In these regions, therefore, a fully 
nonlinear theory must be developed. 

Infinite Ionization Rate Limit 


In the case of infinite ionization rate the 
electron number density is related to the electron 
temperature by Saha's equation (16). Taking the 
square root of equation (16) and expanding the 
right hand side in a power series in T e /\ T e^ we 
obtain for the ratio of the fluctuating density to 
its average value 
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( 55 ) 
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The important thing to be noted is that for de- 
vices of interest in this paper the average elec- 
tron temperature is in general much less than the 
temperature corresponding to the ionization 
energy, i.e. , 


Ce) 
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so that 
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and we have ignored the electron number density 
dependence of the logarithmic factor in the 
Coulomb cross section. 

Expressions for the average electron number 
density and the average of the square of the fluc- 
tuating electron number density can be obtained by 
expanding the electron energy equation in a power 
series in fie/^n e ^. This expansion is carried out 
to third order^in n e /^n e '), neglecting terms of 
the order of T e / Qe) with respect to n e //h e y . 
Since the method is analogous to that considered 
in the orevious limiting case to determine ^T e ^ 
and only the result is presented here. For 

the amplitude function we obtain 
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where 
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Therefore throughout the subsequent analysis terms 
the order of T e ^ T e ^ are neglected with respect 
to fi e /<n e >. 

The equations governing the average current 
density and the fluctuating current density ^ 
are equations (23) and (21) , respectively. These 
equations are specialized to the present limiting 
case by specifying the dependence of d and 6 
upon .T /^T e } and n e /<n e ) in this limit. Ignor- 
ing T e //T e \ with respect to S e /(n e )we obtain 
from equations (6) , (7) , and (8) 
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from equations (19) , (7) , and (8) 
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The coefficient $ 2 , of the linear term is propor- 
tional to the stability term obtained from linear 
wave analysis. Its consequences are reviewed in 
reference (4). It is found that the oscillation is 
most unstable for propagation at the angle given by 


T m 


< 8 > 


VM 2 - 


(63) 


The linear term is then found to be positive and 
hence growing for a Hall parameter 


6 > 6 = 27(1 + K) (A + K) 

c 


The coefficient 


(64) 

is the result obtained 
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in reference (9) and reviewed in reference (4) . 
This result was obtained for an expansion of the 
energy equation to second order in ii e /(n e ^ and 
under the assumption that Z was large in the 
region of interest so that nonlinear. terms not in- 
volving Z could be ignored. These previously 
ignored terms as well as third order terms are 
contained in the coefficient $ 3 . 


tenn 3^, requires large values of (the order 

of the ratio of the mean circumference of the an- 
nular region to its width) for instability to occur. 
Since the analysis indicated that the values of Bq 
in the throat area, where the instability is as- 
sumed to occur, were of order one it does not ap- 
pear that this mode of instability is of importance 
in these devices. 


In figures (7) and ( 8 ) the values of 
a e ff/<°> and 8 e ff/^ 6 ) are plotted, respectively, 
as a function of < 8 > for the theory of refer- 
ence (5), the present third order theory, and for 
the experimental data of reference ( 10 ) for 
(j) = 2 amps/cm2. The theoretical curves are cal- 
culated by the method of reference (9) in which a 
value of the critical Hall parameter, 6 C , is 
chosen to best fit the experimental data of fig- 
ures 7 and 8 . This then allows the energy loss 
factor K to be calculated from equation (63) . 

All other parameters are calculated from the ex- 
perimental data of reference (10). The amplitude 
function, (n|}//n \^, is then calculated from 
equation (56) for 3/3t v* 2 // ( n e ) = 0. 6 e ff/^^ 

and a„£f / (a) are then obtained from equa- 
tions (27) and (28) , respectively, using equa- 
tions (53) and (54). 


It is seen from figures 7 and 8 that the 
additional terms derived in this paper do not have 
a major effect upon the interpretation of the ex- 
perimental data, at least as it is presently in- 
terpreted. The agreement with effective Hall 
parameter is slightly better and the agreement 
with effective conductivity is slightly worse. 

The interesting point, however, is that the terms 
contained in $3 can contribute a negative value 
to the coefficient of the nonlinear term. This 
opens the possibility of instability occurring in 
a region where linear theory predicts stability, 
i.e., 0^ < 0. In this case, instability can occur 
for $3 > $2 the initial amplitude of the 
fluctuation is 




In figure 9 we have plotted (for the conditions of 

, 6 C = 2, Am = A = 1) the 
r an unstable solution in 
0 <_ B c = 2. Obviously 
for these conditions the initial amplitude must be 
large in order for the system to be unstable. One 
means by which an initially large amplitude fluc- 
tuation could occur would be propagation of fluc- 
tuations into the MHD generator from a preionizer. 
This may also be an explanation of the anomolous 
results obtained in the experiments of refer- 
ence (3) during the application of preionization. 


fig. 4 of ref. (5), i.e. 
values of <^>/<n e ) 2 fo 
the linear stable region 


Concluding Remarks 

The nonlinear effects of fluctuations in the 
electron properties of a partially ionized gas have 
teen considered on the basis of a third order 
quasi-linear theory neglecting mode coupling. The 
two limiting cases of frozen flow and infinite 
ionization rate have been investigated with appli- 
cation to MHD generator and high pressure MHD arc 
thrusters operation. 

For MHD generators oscillations occurring in 
the frozen flow limit are damped for normal operat- 
ing conditions. However, where strong preioniza- 
tion is used prior to entering the MHD duct such 
oscillations can be excited and this may explain 
the anomolous fluctuations observed in the experi- 
ments of reference (3). The effect of oscillations 
in the infinite ionization limit has been previ- 
ously studied on the basis of a second order quasi- 
linear theory. (9) a comparison of these results 
to the third order theory developed in this paper 
show little difference in as far as the interpre- 
tation of present experimental data is possible. 
However, the present theory indicates the possi- 
bility of fluctuations occurring in the region 
which is linearly stable provided fluctuations of a 
finite initial amplitude are present in the system. 
In the MHD generator case these may occur as a re- 
sult of strong preionization. 

In high pressure MPD arcs geometrically con- 
fined between two concentric cylindrical electrodes 
it is argued that the most unstable oscillation 
occurs for nearly pure rotation propagation. In 
this case oscillations arising in the infinite 
ionization limit are damped at least for high pres- 
sures. In the frozen flow limit oscillation can 
occur. Based on the assumption of zero gas temper- 
ature in the onset region, which was used in refer- 
ence (2) to reduce the existing experimental data 
for comparison with linear theory, the quasi-linear 
theory predicts unstable growth of the fluctuation. 
Therefore, as a consequence of this assumption, a 
"steady state" fluctuation of finite amplitude is 
not predicted by quasi-linear theory. The validity 
of the zero gas temperature assumption requires 
further experimental information or a more complete 
steady state theory. 


Symbols 


It should be emphasized that the prediction 
of instability in the linearly stable region by 
the present theory is only a possibility since it 
occurs as the nonlinear term begins to dominate. 
Higher order terms may provide sufficient damping 
to stabilize this region. This question can only 
be answered by considering the full nonlinear 
problem. 


A ’K 

defined by eq. 
tively 

(55) and 

(62) , respec- 

A ,a 
ea ea 

defined by eq. 

( 8 ) 


B 

magnetic field 

vector 


£ 

unit vector in 

magnetic 

field direction 


In MPD arc thrusters limited by geometry to C ,C defined by eq. (48) and (49), respec- 

approximately rotational propagation, it has been u 1 tively 

shown in reference ( 2 ) that the linear stability 
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electric field vector 


x,y,z components of x, y, z coordinate system 


e 


magnitude of the charge of an electron 0 


zeroth order value 


I 


electric current density vector 


Superscripts : 


k 


Boltzmann constant 


defined by eq. ( 11 ) 


1 

m 

n 

P 

R 

r_ 

T 

T 

e 

T. 

l 

t 

v 


wave vector 
particle mass 
particle number density 
pressure 

defined by eq. (45) 
position vector 

temperature of neutrals and ions 

temperature of electrons 

temperature equivalent of ionization 
potential 

time coordinate 

gas flow velocity vector 


value in coordinate system moving with 
gas 


Brackets : 

(. . . ) average value of quantity within brackets 
References 

1. Kerrebrock, J. L. , "Nonequilibrium Ionization Due 

to Electron Heating: I. Theory," AIAA Jour - 

nal, Vol. 2, No. 6, June 1964, pp. 1072-1080. 

2. Smith, J. M. , "Electrothermal Instability - An 

Explanation of the MPD Arc Thruster Rotating 
Spoke Phenomena," Paper 69-231, Mar. 1969, 
AIAA, New York, N.Y. 

3. Louis, J. F. , "High Hall Coefficient Experiments 

in a Large Disc Generator," Electricity from 
MHD, 1968 , Vol. II, International Atomic 
Energy Agency, Vienna, 1968, pp. 825-849. 


Z defined by eq. (29) 

6 Hall parameter 


4. Solbes, A., "Instabilities in Non-equilibrium 
M.H.D. Plasmas - A Review," Paper 70-40, Jan. 
1970, AIAA, New York, N.Y. 


v 


total electron momentum collision fre- 5 . 
quency 

ionization coefficient 


Lutz, M. A. , "Radiation and Its Effect on the 
Nonequilibrium Properties of a Seeded Plasma, 
AIAA Journal , Vol. 5, No. 8, Aug. 1967, pp. 
1416-1423. 


ij 


v 

m 


v 

r 


a 
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eff effective value 
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(B) APPROXIMATE GEOMETRY USED FOR ANALYSIS. 
Figure 2. - Thruster geometry. 
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Figure 3. - Ratio of 
effective to zeroth 
order electrical 
conductivity as a 
function of ampli- 
tude function for 
various angles of 
propagation. 
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Figure 4. - Value of tan <p at which j3 eff = 3 0 
as a function of (3 0 . 
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Figure 5. - Amplitude funu 
Hall parameter. 
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Figure 6. - Damping and growth regions for 
rotation disturbance. 
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Figure 7. - Ratio of effective to average conductivity versus 
average Hall parameter. 
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Figure 8. - Ratio of effective to average Hall parameter 
versus average Hall parameter. 
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Figure 9. - Amplitude function 
versus average Hall parameter. 



